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Abstract 

The  electronic  structures  of  LiNi02,  LiA102  andLiNio.5  AI0.5O2  are  investigated  using  density-functional  theory  (DFT)  in  the  local-density 
approximation  (LDA).  The  effect  of  lithium  intercalation  and  the  influence  of  aluminium  doping  on  the  structure  and  electrochemical  prop¬ 
erties  of  LiNi02  are  discussed.  An  increase  in  the  open  circuit  voltage  is  observed  with  A1  doping  in  LiNio.5Alo.5O2  compound. 
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1.  Introduction 

Rechargeable  lithium  batteries  have  been  intensively  stud¬ 
ied  over  the  last  and  the  present  decades  and  are  now  avail¬ 
able  as  competitive  energy  source  [1].  Therefore,  further 
improvements  are  anticipated  concerning  the  environmental 
compatibility,  energy  density,  specific  energy,  cycling  behav¬ 
ior,  etc.  To  achieve  some  of  these  goals,  a  great  number  of 
experimental  works  has  been  devoted  to  the  structure  char¬ 
acteristics  and  the  specific  properties  of  intercalation  com¬ 
pounds  [2-5].  However,  there  still  exist  many  unanswered 
questions  on  the  performance  of  these  materials.  Indeed,  ex¬ 
periments  in  this  area  require  an  immense  amount  of  time 
and  manpower  for  synthesis,  characterization  and  cycling.  In 
parallel,  a  significant  volume  of  theoretical  studies  has  been 
devoted  to  the  calculation  of  band  structure  of  intercalation 
compounds  using  different  computational  methods  [6-10]. 

Lithium  metal  oxide  with  the  layer  structure  LiMCE, 
(M  =  3d  transition  metal)  [11-13]  has  drawn  considerable 
attention  for  lithium-ion  cell  application  due  to  its  high 
energy  density  and  capacity.  It  has,  however,  a  few  disavan- 
tages,  e.g.  stoichiometric  LiNiCE  is  difficult  to  synthesize. 
When  used  as  cathode  material,  LiNiCE  is  less  stable  struc¬ 
turally  upon  cycling  and  shows  severe  capacity  fading,  etc. 
[14,15].  Various  attempts  have  been  made  to  minimize  these 
problems.  Recently,  lithium  metal  oxides  with  a  mixture 
of  cations  on  the  metal  site  (LiMi_>7M/v02)  are  currently 
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of  much  interest  as  rechargeable  battery  cathodes  as  they 
may  offer  the  possibility  to  improve  properties  over  those 
of  single-metal  oxides  (LiM02)  [15-20]. 

Compared  to  the  transition  metals,  A1  has  not  received 
wide  attention  as  a  dopant,  with  only  few  report  on 
LiAlyM  \-y02  (M  =  Co,  Ni)  system  [21-23].  A1  substitu¬ 
tion  of  the  transition-metal  cation  has  been  shown  theoret¬ 
ically  [6-10]  and  experimentally  [2-5]  to  increase  the  cell 
voltage.  This  effect  has  been  verified  for  LiAlvCoi_>702 
solid  solution  [24].  Moreover,  LiAlCE  is  stable  in  the 
Q'-NaFe02  structure  at  temperature  below  ~600°C  [25] 
and  it  has  a  stabilizing  effect  of  the  structure  as  reported  for 
monoclinic  LiAlyMni-^CE  [26].  Its  low  cost  and  low  den¬ 
sity  makes  LiA102  attractive  as  an  intercalation  compounds 
constituent. 

The  important  features  of  a  lithium/metal-oxide  com¬ 
pound  for  battery  applications  include  the  voltage  at  which 
it  exchanges  lithium,  the  amount  of  lithium  that  can  be 
reversibly  intercalated  and  the  stability  of  the  material.  The 
first  two  properties  determine  the  energy  density  and  the 
latter  limits  the  lifetime  of  the  battery.  Ceder  et  al.  [24,27] 
have  demonstrated  that  the  average  potential  for  intercala¬ 
tion  can  be  obtained  from  the  energy  change  in  the  reaction: 

Li  +  MxM[_x02  ->  UMxM[_x02  (1) 

First-principles  computation  is  able  to  calculate  the  energy 
of  the  three  compounds  in  Eq.  (1),  and  hence  the  average 
potential,  for  any  metal  M  or  combination  of  them,  in  any 
structure,  whether  these  conditions  have  been  experimentally 
realized  or  not. 
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In  this  paper,  the  Al  effect  on  the  structure  and  elec¬ 
trochemical  properties  of  LiNiC>2  compound  has  been  in¬ 
vestigated.  Numerical  calculations  were  performed  on  the 
LiNiC>2,  LiA102  and  extended  to  quaternary  spinel  phase 
LiAlo.5Nio.5O2. 

2.  Computational  details 

Our  calculations  were  performed  using  density-functional 
theory  (DFT)  [28-30]  with  the  exchange-correlation  func¬ 
tional  treated  in  the  local-density  approximation  (LDA) 
as  parametrized  by  Hedin  and  Lundqvist  [31,32].  The 
eigenvalue  problem  was  solved  applying  the  full-potential 
linearized  augmented-plane-wave  (FP-LAPW)  method  as 
implemented  in  the  WIEN2K  code  [33-35].  This  method 
has  been  known  to  be  one  of  the  most  accurate  quantum 
mechanical  methods  since  it  improves  upon  linearization 
and  makes  possible  a  consistent  treatment  of  semicore  and 
valence  in  one  energy  windows  and  allows  inclusion  of  local 
orbits  in  basis.  In  order  to  achieve  a  satisfactory  degree  of 
convergence,  the  cutoff  parameters  used  in  the  LAPW  calcu¬ 
lation  are  summarized  shortly  in  the  following:  in  the  muffin 
tin  spheres  the  wave  functions  were  expanded  in  spherical 
harmonics  with  an  angular  momentum  up  to  /^x  =  10.  Non- 
spherical  contributions  to  the  electron  density  and  potential 
within  the  MTs  were  considered  up  to  /^ax  =  4.  The  cutoff 
for  the  Fourier- series  expansion  of  the  interstitial  electron 
density  and  potential  was  chosen  to  be  Gmax  =  14bohr_1 
and  the  cutoff  for  the  plane  wave  basis  set  RKmSLX  =  8.  In 
general,  the  local-density  approximation  presents  some  de¬ 
ficiencies  as  a  modest  increase  in  computational  workload. 

To  calculate  the  electronic  structure  of  the  quaternary 
spinel  phase,  a  layered  trigonal  supercell  was  used  which 
contains  one  nickel  and  aluminium  ions,  two  lithium  and 
four  oxygen  ions,  thereby  providing  the  LiAlo.5Nio.5O2  con¬ 
tent.  In  such  compound,  only  three  points  are  used  to  rep¬ 
resent  the  Brillouin  zone  comparatively  to  ten  points  for 
LiNi02  and  LiA102. 

To  determine  the  average  intercalation  voltage  (AIV), 
the  total  energy  of  three  structures  is  required  using  the 
FP-LAPW  method:  metallic  lithium  in  body  centered  cu¬ 
bic  and  the  oxide  structures  before  and  after  intercalation. 


The  energies  of  these  oxides  are  calculated  in  the  a-NaFe02 
structure.  The  symmetry  of  this  structure  is  R3m.  Three  de¬ 
grees  of  freedom  exist:  the  a-  and  c- lattice  parameters  and 
the  internal  coordinate  z.  Lattice  constants  and  internal  co¬ 
ordinate  are  relaxed  to  find  the  lowest  energy.  Combining 
the  basic  thermodynamics  [36-38]  and  these  total  energies, 
the  AIV  of  these  compounds  can  be  predicted. 

For  the  electronic  configuration,  the  states  treated  as  bands 
were:  Li(2s),  Ni  (4s,  3p,  3d  ),  Al  (3s,  3p)  and  O  (2s,  2p). 
The  muffin-tin  radii  chosen  are  1.8,  1.9,  1.9  and  1.55  a.u.  for 
Li,  Ni,  Al  and  O,  respectively. 

3.  Result  and  discussions 

3.1.  Structural  properties 

Using  the  first-principles  all  electron  full  potential  LAPW 
method,  the  geometry  optimizations  were  performed  on  the 
layered  trigonal  LiNi02,  LiA102  and  LiNio.5Alo.5O2  as  well 
as  the  corresponding  host  structure.  The  results  for  the  a  and 
c  lattice  constants,  the  oxygen  positional  parameter  z  and  the 
volume  are  listed  in  Table  1 .  Where  available,  experimental 
data  is  added  in  parenthesis. 

For  LiNi02,  the  a  and  c  lattice  constants  are  approxi¬ 
mately  2%  smaller  than  Ni02.  In  the  case  of  the  LiA102,  the 
a  and  c  lattice  constants  are  respectively  1  and  6%  greater 
than  the  delithiated  structure.  The  effect  of  the  aluminum 
doping  on  the  lithium  nickel  oxide  has  been  investigated 
by  substituting  50%  of  the  nickel  by  aluminium  atoms.  The 
a-parameter  is  slightly  increased  by  about  1%  and  an  im¬ 
portant  increase  (15%)  of  the  c-parameter  is  observed  when 
lithium  is  intercalated  into  the  host  structure.  As  is  typical 
in  LDA,  the  optimized  volumes  from  our  calculations  are 
smaller  than  the  experimental  values. 

The  a  lattice  constant,  in  all  intercalated  compounds,  is 
slightly  modified  by  the  lithium  intercalation  and  this  vari¬ 
ation  can  be  directly  related  to  the  metal-metal  distance.  In 
the  Table  2,  the  metal  to  metal  distance  and  lithium-oxygen 
bond  lengths  for  the  intercalated  and  deintercalated  com¬ 
pounds  as  well  as  the  metal-oxygen  bond  lengths  are 
given.  The  hexagonal  layered  structure  can  be  described 
by  a  slab  formed  by  a  plane  of  transition  metal  atoms,  in 


Table  1 


Computed  crystallographic  (the  lattice  constants  a  and  c,  and  the  oxygen  positional  parameter  z)  parameters  for  LLN1O2,  LLAIO2,  and  LLNio.5Alo.5O2 
compounds  in  the  a-NaFe02  structure  (with  2c  =  0  and  1) 


Compounds 

a  (A) 

c  (A) 

z 

v  (A3) 

□Ni02 

2.79,  2.87  (2.84) 

13.13,  11.70  (13.47) 

0.263,  0.255 

29.61,  27.82  (31.38) 

LiNi02 

2.82,  2.93  (2.88) 

14.06,  13.20  (14.19) 

0.259,  0.260  (0.259) 

32.30,  32.71  (33.97) 

□aic>2 

2.85  (2.87) 

11.33  (11.70) 

0.253  (0.255) 

26.69  (27.82) 

LiA102 

2.78  (2.8) 

13.32  (14.23) 

0.258 

29.70  (32.21) 

DNio.5Alo.5O2 

2.98 

12.58 

0.257 

32.28 

LiNio.5Alo.5O2 

3.01  (2.855) 

14.84  (14.24) 

0.254 

38.95  (33.62) 

The  □  indicates  the  host  structure  without  intercalation.  The  □  correspond  to  2c  =  0. 
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Table  2 


Computed  average  intercalation  voltage,  in  volts,  for  LiMC>2  com¬ 
pounds,  and  metal-oxygen,  sodium-oxygen  and  metal-metal  distances  in 
(□/Li)MC>2  compounds,  all  entries  in  (A) 


LiNi02 

LiA102 

LiNio.5  AI0.5O2 

AIV’s 

3.32  (3.9  [11]) 

5.3  (5.4  [24]) 

4.0 

rf(M-O)  (D/Li) 

1.85/1.93 

1.88/1.89 

1.62/1.68 

d(U-0)  (D/Li) 

-/2.08 

-/2.01 

—/1 .70 

d(M-M)  (D/Li) 

2.79/2.82 

2.85/2.78 

2.98/3.01 

The  □  indicates  the  host  structure  without  intercalation.  The  □  correspond 
to  v  =  0. 


between  two  planes  of  oxygen.  Interaction  between  atoms 
in  a  slab  (M-O)  are  mainly  covalent,  while  the  slab- to- slab 
interactions  (Li-O)  are  generally  weaker.  The  M-0  dis¬ 
tance  (shorter  than  the  Li-0  distance)  decreases  from  Ni 
to  Nio.5Alo.5-  The  covalent  character  appears  to  be  more 
important  in  the  Nio.5Alo.5  in  comparison  to  the  Al  and  Ni 
compounds.  This  theoretical  result  confirms  the  important 
role  of  the  aluminium  effect  in  order  to  stabilize  the  lithium 
nickel  oxide  in  the  a-NaFe02  structure  and  corroborates 
experimental  observations  [2-5]. 

When  Li  is  intercalated,  the  M-0  bond  length  increases 
in  Ni  and  Nio.5  AI0.5  compounds.  This  structural  effect  is  due 
to  the  lithium  influence  into  the  space  between  the  O-M-O 
slabs,  which  pulls  the  oxygen  towards  it  and  leading  an 
increase  of  the  M-0  distance.  For  Al  compound,  same  M-0 
bond  length  is  reported.  This  feature  is  attributed  to  the  small 
nuclear  Al  charge  which  induces  an  effectively  smaller  metal 
ion. 

Table  2  gives  the  predicted  average  intercalation  voltage 
for  Ni,  Al  and  Nio.5Alo.5  compounds.  For  lithium  insertion 
in  spinel  structures,  it  can  be  seen  that  there  is  a  reason¬ 
able  agreement  between  experimental  and  calculated  values, 
with  the  experimentally  potential  consistently  being  under¬ 
estimated  by  about  by  about  1/2  V  for  LiNi02  compound. 
For  LiA102,  the  AIV  is  in  the  same  order  with  the  theoreti¬ 
cal  value  reported  by  Aydinol  et  al.  [38].  Moreover,  the  Al 
compound  induces  clearly  an  AIV  which  is  higher  than  that 
of  any  lithium  transition-metal  oxide.  From  the  value  of  the 
average  intercalation  voltage  obtained  for  the  mixed  com¬ 
pound  (LiNio.5  AI0.5O2),  it  is  clearly  that  the  potential  in¬ 
creases  systematically  with  the  substitution  of  the  nickel  by 
the  aluminium.  This  variation  confirmed  the  positive  effect 
of  aluminium  as  it  is  demonstrated  experimentally  [2-5]  and 
predicted  theoretically  [6-10].  In  this  case,  the  AIV  increases 
by  about  0.7  V  in  comparison  to  pur  nickel  compound. 

3.2.  Electronic  properties 
3.2.1.  Band  structure 

In  Fig.  la  and  b,  the  LDA  calculations  of  the  band  struc¬ 
tures  of  LiNiC>2  and  LiA102  as  well  as  the  delithiated  corre¬ 
sponding  compounds  in  the  a-NaFe02  structure  are  shown. 
The  band  structures  and  band  ordering  are  very  similar  in 
the  first  region  characterized  by  an  energetically  low  lying 


bands  below  the  Fermi  level.  These  bands,  derived  from  the 
2s  states  of  the  oxygen  atoms,  do  not  contribute  to  the  bond¬ 
ing.  For  the  aluminium  compounds,  the  second  region  is 
composed  of  six  bands  mainly  due  to  2p-0.  In  the  case  of 
the  nickel  compounds,  nine  bands  are  obtained  and  they  are 
originated  mainly  from  2p-0  in  some  part  of  the  Brillouin 
zone  and  3d-Ni  in  another  part  of  the  zone,  as  can  be  seen  by 
an  overlap  and  a  mixing  between  these  bands  in  the  T -point 
(six  bands  are  the  2p-0  character  and  three  bands  are  the 
d(t2g)-Ni  character).  The  band  ordering  is  really  different  at 
the  T-point  in  NiC>2  comparatively  to  LiNi02-  In  Ni02,  the 
first  three  bands  are  originated  from  2p-0  with  an  important 
contribution  of  the  d(^g)-Ni  which  implies  the  presence  of 
strong  a  bond  without  excluding  the  n  bond.  The  second 
three  bands  are  dominated  by  the  d(^g)-Ni  bands  and  the 
last  are  derived  from  2p  states  of  the  oxygen.  For  LiNiC>2, 
the  first  part  remains  the  same  with  an  increasing  of  the  par¬ 
tial  charges  d(eg)-Ni  and  a  decreasing  of  d(t2g)-Ni  which  in¬ 
dicates  the  presence  of  very  strong  a  and  reduced  n  bonds. 
These  bands  are  followed  by  p-0  bands.  In  the  third  order, 
the  d(t2g)-Ni  bands  are  obtained.  The  highest  band  shown  in 
Fig.  la,  are  originated  from  d(eg)-Ni  in  the  NiC>2  and  LiNiC>2 
compounds.  From  Fig.  lb,  the  highest  band  above  Fermi 
level  is  originated  from  s-Al,  because  in  AIO2  and  LiA102, 
the  repulsive  interaction  with  the  O  s-bands  shifts  this  band 
to  higher  energies  above  the  E p.  This  band  is  followed  by 
p-Al  band  which  does  not  contribute  to  the  bonding. 

Comparing  the  band  structure  of  the  intercalated  and 
deintercalated  compounds,  all  the  bands  are  shifted  to  higher 
energies  due  to  the  lithium  intercalation.  The  Fermi  energy 
increases  with  the  increasing  of  the  number  of  valence 
electron  from  deintercalated  to  intercalated  compound.  Fi¬ 
nally,  the  intercalated  compound  is  more  covalent  due  to 
the  increasing  of  the  o  bond  comparatively  to  the  parent 
compound. 

The  electronic  structure  of  the  quaternary  alloy  and  the 
corresponding  host  structure  is  similar  to  those  obtained  for 
the  oxide  compounds.  The  narrow  band  which  has  little  in¬ 
teraction  with  other  states  is  dominated  by  the  s-0  charac¬ 
ter.  The  next  bands  consist  of  2p-0  and  lower  3d(/L2g)-Ni 
orbitals.  In  these  bands,  the  Ni-0  mixing  is  also  shown.  The 
band  at  the  Fermi  level  is  mainly  3d(^g)-Ni  character  fol¬ 
lowing  by  the  s-Al  and  p-Al,  respectively.  As  Li  is  inter¬ 
calated  into  the  host,  the  Fermi  level  moves  up  in  energy 
(3.54-4.80  eV)  and  the  all  bands  are  shifted  up.  Moreover, 
the  band  ordering  is  similar  comparatively  to  the  host  struc¬ 
ture. 

3.2.2.  The  density  of  states 

The  total  DOS  are  shown  in  Fig.  2  for  LivNixAli_x02, 
where  v  =  0;  1  and  y  =  0;  1.  The  lithium  intercalation  effect 
can  be  examined  into  nickel  compound  by  comparing  the 
total  and  partial  DOS  for  LiNi02  and  Ni02  shown  in  Fig.  2a. 
It  is  clear  that  the  0-2s  band  dominates  the  first  region  in 
the  lithiated  and  delithiated  compound.  In  the  second  region 
for  the  LiNi02,  the  0-2p  band  extends  over  the  wide  energy 
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Fig.  1.  The  LDA  energy  band  structure  of  LiNixAli_x02  compounds  in  the  layered  structure,  where  x  =  0,  1  and  0.5. 


range  with  a  important  contribution  of  the  Ni-3d(£g)  bands.  0-2p  and  Ni-3d  bands  in  NiC>2  has  reduced  to  a  minimum 

Moreover,  the  Ni-3d  band  is  mainly  located  in  the  higher  of  the  density  of  states  which  separe  these  two  bands  in 

energy  region.  The  lithium  intercalation  into  NiC>2  makes  LiNi02.  In  addition,  the  Fermi  level  Ep  shift  up,  from  the  top 

the  0-2p  band  to  shift  up  largely  (from  —0.24  to  —0.10  Ryd  of  the  dfeg)  to  the  bottom  of  the  d(eg)  band  in  the  LiNiC>2. 

in  Ni02  and  LiNi02,  respectively).  The  overlap  between  the  This  feature  is  attributed  to  the  fact  that  the  intercalated 
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electron  has  to  be  accommodated  in  the  high  energy  Ni-d(^g) 
band.  For  LiA102  and  AIO2,  (Fig.  2b),  the  band  ordering  is 
very  similar  for  the  first  and  the  second  region.  The  main 
difference  concerns  the  third  region  which  is  constituted  to 
the  Al-3s  and  Al-3p  characters. 

3.2.3.  Valence  electron  density 

To  study  the  Li  intercalation  effect  on  the  valence 
electron  density,  we  have  calculated  the  difference  be¬ 
tween  the  valence  electron  density  in  LiNi02,  LiA102  and 


LiNio.5Alo.5O2.  For  this  purpose,  we  have  chosen  the  (110) 
crystallographic  plane  which  contains  lithium,  oxygen  and 
nickel  (aluminium)  ions  for  all  compounds.  For  these  cal¬ 
culations,  all  atomic  positions  were  assumed  to  remain  the 
same  after  lithiation. 

From  Fig.  3a,  it  seems  that  there  is  no  doubt  that  the 
charge  transfer  is  from  the  Li  to  the  Ni  and  O  ions  in 
LiNi02  and  only  to  O  ions  in  LiAlC>2.  For  the  quaternary 
alloy  (LiNio.5Alo.5O2),  the  same  feature  is  observed.  From 
Table  3,  the  difference  in  the  electronic  charge  inside  the 
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(C) 


Fig.  3.  Difference  in  the  charge  density  between  the  lithiated  and  the 
unlithiated  materials.  The  plane  shown  is  parallel  to  the  c-axis  and  con¬ 
tains  lithium,  oxygen  and  nickel  (aliminium)  ions,  (a)  LiNi02/Ni02,  (b) 
LiA102/A102,  and  (c)  LiNio.5Alo.5O2/Nio.5Alo.5O2. 


Table  3 


Electron  charge-transfer  values  to  oxygen  (O)  and  transition  metals  (M) 
upon  intercalation  of  Li  into  MO2  compounds  in  the  a-NaFe02  structure 


Ni-O 

Al-O 

(Nio.sAlo.5)-0 

M 

0.21 

— 

0.31  (Ni) 

-(Al) 

O 

0.21 

0.36 

0.154  (Ol) 

0.227  (02) 

Note  that  there  are  two  anions  per  unit  cell. 


atomic  sphere  between  the  crystal  (LAPW)  value  G^iystd 

LiNixAli_x02  and  the  superposed  atomic  value  Q super6  in 
NixAli_*02  (where  x  =  0,  1  and  0.5),  which  is  the  quantity 
defined  the  charge  transfer  in  a  compounds  before  and  after 
the  intercalation  is  presented.  From  these  values,  a  charge 
transfer  of  about  0.2  an  electron  from  the  lithium  to  the  nickel 
and  the  oxygen  atoms  occurs  in  the  Nickel  compound.  This 
result  is  similar  to  those  obtained  by  Aydinol  et  al.  [37,38]. 
In  the  case  of  the  aluminium,  there  is  exclusively  a  charge 
transfer  to  the  oxygen  atoms  of  about  0.36  electron.  For  the 
quaternary  alloy,  the  increasing  of  the  charge  transfer  to  the 
nickel  by  about  32.5%  is  observed  comparatively  to  LiNi02. 
Moreover,  the  charge  transfer  to  oxygen  is  slightly  increased 
in  comparison  to  pur  compound  (LiNiC>2  and  LiA102). 

4.  Conclusion 

In  this  paper,  we  showed  that  the  linearized  augmented 
plane  waves  method  can  be  used  to  estimate  the  average  in¬ 
sertion  potentials  as  well  as  the  electronic  structure  of  Ni02, 
AIO2  and  Nio.5Alo.5O2  cathodes  for  secondary  lithium  bat¬ 
teries.  The  agreement  between  the  experimentally  mesured 
and  the  calculated  average  potential  is  consistently  underes¬ 
timated  by  about  1/2  V  for  LiNi02,  whith  confirms  the  ac¬ 
curate  prediction  of  the  LAPW  method  for  all  compounds. 
Moreover,  we  find  that  the  substitution  of  Ni  by  Al  lead 
to  higher  Li  intercalation  voltage  comparatively  to  the  pur 
compound  (LiNi02)  and  stabilize  the  structure.  We  have 
shown  that  the  charge  transfer  from  Li  to  Ni  atoms  has  an 
effect  on  the  AIV  without  exclude  the  great  importance  of 
the  charge  transfer  from  Li  to  O  atoms  as  demonstrated  by 
many  authors. 
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